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This technical paper describes an analysis tool that utilizes a database of transfer 
functions (TFs) to generate quick turn-around vibroaeoustic response estimates suitable for 
developing initial vibration environment vendor specifications during the early stages of a 
new launch vehicle design program. The tool draws TFs from a database, combines TFs, and 
multiplies these by input excitations to estimate vibration responses. Initially, the database is 
populated with two sets of uncoupled TFs, the first set representing the vehicle panel 
unloaded vibration response TF, designated as H“, and the second set representing the free- 
free component equipment (component only, no vehicle panel) vibration response TF, 
designated as H p . For a particular configuration undergoing analysis, the appropriate H“ 
and H p , are selected and coupled to generate an integrated TF, designated as H“ +p . This 
integrated TF, H u+|! , is then used with the appropriate input excitations to estimate vibration 
responses. This simple, but powerful, tool enables a user to estimate vibration responses 
without directly using finite element models, as long as suitable H“ and H p sets are defined in 
the database libraries. The paper will discuss the preparation of databases and provide the 
assumptions and methodologies necessary to combine H“ and H p sets into an integrated H u+P . 
Validation of the coupled TFs will also be presented. 
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NESC = NASA Engineering and Safety Center 
RMM = Response Matching Method 
SPL = Sound Pressure Level 

SPCD = Secondary Component Loads Development Program 
Transfer Function Designations: 

H“, - Exterior Vehicle Panel Unloaded Response Transfer Function defined at Interface Response Locations 
- Free-Free Component Equipment Response Transfer Function from Interface Base Excitations 
H" +,i - Integrated System Response Transfer Function from excitation over surface of Exterior Vehicle Panel 
II /pi ( co ) is the transfer function between acceleration at point j and pressure p h on patch b , 

C, n is the critical damping ratio for mode m , 

F. is the static force at point k associated within a unit pressure on patch b , 
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M J - Mode shapes for modes 1 M for bare skin at a few reference locations j , 
M ^ - Natural frequencies for modes 1 — ^ of bare skin, 

- - Mode shapes for modes 1 ~ M fo r component-loaded skin at locations ^ , 

- Modal forces for modes 1 — M f or component-loaded skin at locations ^ , 


¥ 


qM 




Natural frequencies for modes 1 M 0 f component-loaded skin. 


<j) jm is the ni mass-normalized mode shape at response point j , 

(f) km is the mass-normalized mode shape at point k in the pressure patch, 
co is the circular frequency, 

C0 m is the circular natural frequency of mode rn , 

M is the number of retained modes. 


Static Data: 

F b=l F bi F bi ••• F bN b \ , 

b = l,2,-N 0 


Force distribution on GRIDs for unit pressure on patches ^ , 


^ bc X,:: ^ + ^ + ^ Zb Z<r ^ - Distance between CGs of patches ^ and c . 


Pressure Data: 


P f (of) 

re] - Reference pressure autospectrum. 

N b is the number of GRIDs in the pressure patch. 
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W ( co ) 

bbK - Scaling functions for non-uniform pressure autospectra over the entire skin. Calculate these to coincide 

W 

with the zone autospectra specified in the Vibroacoustic Loads Databook, if available. That is, use the same bb on 

W 

patches within specified launch vehicle zones. Set bb to unity if unknown. 


I. Introduction 

P ROVIDING suitable vibration requirements for launch vehicle equipment is a challenge faced by each new 
program. NASA’s Engineering and Safety Center (NESC) has identified methodology for producing vibration 
response estimates for mass loaded vehicle panels as one of the top risks for new vehicle programs, observing that 
several different methodologies are in use that may generate differing response predictions. Program managers need 
realistic estimates for vibroacoustic response environments so that vibration environment specifications can be 
provided for the design of hardware. Accurate prediction of vibration environments in early stages of a new 
program helps to assure adequacy of the hardware to withstand the expected environments without specifying 
overly-severe environments. Specifying environments that are overly-conservative can drive up costs and add 
weight to the vehicle design. Providing realistic vibration environments can reduce schedule impacts, as the 
hardware does not require extreme design measures to withstand the vibration environments and environmental tests 
may be easier to perform. Imposing less severe test requirements may expand the number of test facilities that are 
able to satisfy required test levels. Engineering organizations engaged in a new vehicle design need a process that is 
quick to implement and easy to use. A complicated approach may present a schedule that it would be wiser to 
avoid. Therefore, engineering would like a proven process that provides seamless cross-cutting capability. The aim 
of this paper is to make such a methodology more accessible, as a capability, to the entire vibroacoustic loads 
community. Providing a tool that produces reliable vibration response estimates and is easy to use is the motivation 
for this work. 


II. Background 

In response to the NESC’s identification of a critical technological need, the authors have presented validated 
methodologies for calculating both vibration response and dynamic loads for equipment mounted to vehicle exterior 
panels. 9 The validation was based on acoustic response tests of a flight-like vehicle panel configured with different 
equipment masses. 10 ' 11 ' 12 We are seeking to make that demonstrated methodology more accessible as a capability to 
the wider loads community. 


III. Innovations for the End User 

The motivation for development of this tool is to address both the need and the vision to provide “vibroacoustic 
response estimates,” using the most accurate game-changing estimation techniques, while minimizing complexities 
for the end-user. While a more straight-forward, but limiting, approach could be utilized that uses just a single 
transfer function representing a specific case, the goal of the proposed versatile approach is to provide a multi- 
transfer function capability that enables an end-user to address and investigate several trade configurations with 
varying complexities. 

The proposed approach is intended to be simple for the end-user, no longer requiring the end-user to work with 
finite element models. This capability is accomplished by leveraging off the work of others that have previously 
prepared the models and entered the response characteristics of these models into the database library. Since the 
unloaded panel and the component equipment transfer functions are generated separately, the end-user has numerous 
integrated mass-loaded panel options available. Integrated transfer functions are generated automatically in the tool 
by coupling the unloaded panel and component equipment panel transfer functions using the numerical procedure 
outlined herein. The end-user makes the following menu selections: 


9 Peck, J., p 5-21. 

10 Smith, A.,. 

11 Driskill, T..,. 

12 Rodgers, C.,. 
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1 . The appropriate H“ set - which depends on: 

a. Zonal Exterior Vehicle Panel 

b. Number and location of Interfacing connections to Single Equipment Assembly 

2. The appropriate H |! set 

a. Number and location of Interfacing connections of Single Equipment Assembly to primary vehicle 
panel 

b. Equipment Assembly Mass and CG. 

c. Estimated Damping Spectrum 

3. The appropriate SPL description of the Excitation Environment. 
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Figure 1. Process flow for Database Tool Preparation Team and Vibration Environment Preparation Team 


From these simple inputs, the Database Tool constructs the H“ +|3 TFs for the integrated system at each interface 
location. Next, narrowband vibration PSDs are calculated using the appropriate Pressure Auto- and Cross-Spectral 
densities that correspond to a selected zonal SPL at the interface nodes and provides them as a narrow band 
spectrum. 

In order to provide the Max Expected Environment, the end-user may only need to construct a smoothed 
envelope of the narrow band PSDs output from the analysis tool, adding an appropriate uncertainty factor to result in 
vibration requirements suitable for a vendor specification. 
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The calculated environment is intended to provide greater accuracy because more comprehensive design 
details are represented in the calculations than are typically considered using first-order estimates developed 
using traditional empirical scaling and mass attenuation techniques. 


Figure 1 presents the expected process flow diagram to develop the tool, recommending two teams work in 
parallel. Efficient use of engineering resources may be achieved by assigning a small group to populate the 
Database Tool, “The Preparation Team,” while assigning others to work on the tool itself, incorporating 
formulations, generating a GUI, and generating enveloping routines for requirement specification, “The Production 
Team.” Since detailed component equipment interface details are not available until long after vendor specifications 
are released for a new vehicle program, development of a generic set of H is advised. These can be developed long 
before the actual component installation design details are decided. The generic transfer function set should include 
options that are realistic regarding the number and spacing of component equipment interfaces with the primary 
vehicle structure exterior panels. Preparing H 1 ' for a wide range of weights and CG offsets is recommended. Once 
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Figure 2. Notes on Setting up Uncoupled Transfer Function Ca 


the primary structure design begins to solidify, a set of FT TF would be prepared using the Acceleration/Pressure 
Transfer Function (APTF) method with a unit P le f for the first Vehicle Zone (see equations 10 and 11). Guidance on 
calculating uncoupled TF is presented in Figure 2. 

Equipped with generic FI |! and vehicle design specific 1-1“ sets, the Production Team would begin to work in 
parallel to the preparation team. They would be able to provide both: 

1. The Zonal Vibration Environments necessary as inputs to develop the Vehicle Program-Loads and 
Environments Data Book 
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2. The vibration requirements necessary to issue vendor specifications for 
component design. 

Figure 3 provides guidelines for boundary conditions to use in calculating the 
vehicle panel side transfer functions. The vibration environment estimates should be 
more accurate than those developed using the heritage approach, especially in 
important frequency ranges. A reduction in the weight and complexity of 
components and their attachment brackets may be possible by producing vibration 
environment estimates with greater accuracy. 


IV. Basis Methodologies 

A pair of methodologies was previously presented by the Authors. 1 ’ The first 
method is the APTF, which applies a correlated pressure field across the surface of a 
vehicle panel. In this application, a Diffuse Acoustic Field (DAF) pressure 
assumption is adopted. Pressure spectra provided as “RMS sound pressure levels” 
are first converted to “pressure autospectral density”. Then cross-spectra associated 
with the pressure field excitation are calculated according to the best fit of for a DAF 
using Equation (4). To implement this method, the excited vehicle panel surface 
must be subdivided into pressure patches. The distance between the center of each 
patch is selected based on the wave number required to calculate cross-spectra 
phasing between each patch. 

An acceleration response transfer function is first described in terms of an 
excitation on a single patch in Equation (5). But the transfer function is fully 
developed to include the spatial correlation of the pressure field across the many 
patches of the vehicle panel in eq (10). 

The second method, the Response Matching Method (RMM), may be applied when the vibration response of a 
bare vehicle panel is a known input to the problem. Several exterior vehicle panel response transfer functions, each 
defined at appropriate interface response locations, 14“, can be obtained using a portion of eq (10). Using the finite 
element method it is also possible to calculate a corresponding transfer function, for the integrated system (the 
same panel with equipment integrated on it) from eq (14). 

One of the efforts involved in development of a working database tool will be to demonstrate the process of 
calculating the H al|! , without using an integrated FEM of the exterior vehicle panel and the supported equipment. 
Instead, the H“ |! is produced using the analysis tool, employing component mode synthesis from the union of 
selected H“ and H 1 ' options available within the library database. 

The equations for the validated and very effective APTF and RMM approaches may be understood from the 
following equations and discussion. A random pressure field on a launch vehicle skin surface may be approximated 
by dividing the surface into “patches,” or regions of uniform pressure with no phase offsets within the patch. The 
size of each patch must be chosen to be small enough to justify the assumption of uniform (though dynamically 
varying) pressure with zero phase offset across the patch. The selection of the appropriate patch and element size is 
not included in this monograph, but is addressed in a separate study. 

The pressure may be defined as a stationary Gaussian random field with spatially varying autospectral density. 
The pressures on any pair of patches may be correlated, exhibiting a non-zero cross-spectral density between them. 

The random pressure field is thus a square Hermitian matrix of spectral densities of dimension / ^/> , the total number 
of pressure patches. The pressure autospectra occur on the diagonal of the matrix. The cross-spectra includes the off- 
diagonal terms. The random pressure field on all patches may be written as 


Figure 3. A Useful Boundary 
Condition for Vehicle Panel Might be 
an Elastic Model of Structures 
Forward and Aft. Including a 360 
Cylindrical Section. Using 
Symmetric Boundary Conditions may 
be an Another Option. 
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’ Peck, J., p 16-21. 
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where P hc = P, , and the asterisk denotes the conjugate operator. If spatial functions yip), R) are defined that 
relate the autospectra to the cross-spectra, eq (1) may be written as 
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where P hc = yjP bb P cc and y bh have been added to the diagonals for generalization. The expression for P hc arises 
from an inequality requirement on the coherence, which states that 


0 < 


pjm r 

PmMPJoj) 


< 1.0 


(3) 


For a diffuse field, the spatial fimctions y may be expressed as 


r b c(°p R bc) 


sin(i? fe r(6>)) 

R bc K{(D) 


(4) 


where R hc is the distance between the area CGs of patches b and c , fc(co) = 6)/ C a , and C Q is the speed of 
sound through the fluid medium adjacent to the patch material. The patch CGs are constrained to lie on the curved 
skin surfaces. When b = C , the spatial functions coincide with the patch autospectra, the distance R between 
patches vanishes and y bh — > 1.0 in the limit as R hh — > 0 (L’Hopital’s Rule). Finally, the patch autospectra may 

be expressed as products of frequency-dependent scaling functions W bb io) and an arbitrary reference 
autospectrum, (e.g., one of the patch autospectra could be selected, but this is not required.) 

P bb (co) = W bb (co)P ref (co ) (5) 


Substituting eq (5) into eq (2) results in an expression for the diffuse-field pressure model used for the lift-off 
environments (Cases 01 - 08) in this study: 
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where W bc (co) = ■\j^ T hh (co) W cc ( CO) . It is important to recognize that all of the pressure auto and cross-spectra 
represented in eqs (1) - (6) are known, and are related to the reference spectrum P re/ through the scaling functions. 
If the scaling functions are unknown, then the skin structure under consideration must be limited to regions with 
similar pressure autospectra and W bb must be set to unity. The scaling functions may also be grouped into fewer 

distinct sets for large regions showing uniform pressure autospectra. That is, many of the functions may be assigned 
the same value for various zones of the launch vehicle even if the zone is subdivided into many patches. 

Now consider the acceleration frequency response a u (co) at a point j on the skin to a unit oscillating pressure 
p n (co) = COS (cot) on patch b . The frequency response is also the transfer function between any acceleration 
response Cl ■ ( CO) at j to arbitrary pressure p h (co) on patch b. Note that the response point locations are 

unrestricted; they may be located either within or external to the patch. In NASTRAN, the unit pressure is 
distributed as a set of non-uniform forces on the grid points comprising the shell elements within the patch. This 
non-uniform force distribution for a unit pressure may be obtained with an OLOAD request at the patch GRIDs in 
NASTRAN SOL 101 (linear static loads solution). The force distribution may then be swept through the frequency 
range of interest to obtain the frequency response. The acceleration/pressure transfer function for a single patch may 
be expressed explicitly as the sum of weighted acceleration/force transfer functions. The weighting factor is the 
static force F k at each input location k on the patch due to a unit pressure (from an OLOAD request in SOL 101): 


N b M 


k = 1 m = 1 


~ 0j2 </> Jm Am 


co] - co +i co m co 

m ~ m m 


(V) 


where 


II lpi (co) is the transfer function between acceleration at point j and pressure p h on patch b , 

F k is the static force at point k associated with a unit pressure on patch b , 

/ . th • • • 

(p j m is the ill mass-normalized mode shape at response point j , 

<f> km is the fil' h mass-normalized mode shape at point k in the pressure patch, 
co is the circular frequency, 

CO m is the circular natural frequency of mode m , 

C, m is the critical damping ratio for mode m , 

K is the number of GRIDs in the pressure patch, 

M is the number of retained modes. 


The acceleration PSD response Aj b (co) to a random pressure P bh (co) on patch b is the squared magnitude of 
the acceleration/pressure transfer function in (7) multiplied by the pressure PSD, or 
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The total response at location j includes the autospectra from the pressures on all of the patches and also from 
non-zero cross-spectra between any two patches, or 


b c*b 


4 w = f p » + 
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b c 
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Expressing eq (9) in terms of the reference spectrum and the spatially dependent cross-spectra of eqs (4) and (6), 
we obtain 


_ p _ p _ p _p sinf k R 1 
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Note that the spatial functions ^ reduce to unity for b C ; as mentioned previously. Eq (10) may be expressed 
in matrix form for computational efficiency in Matlab as 
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Where ' denotes the Hermitian conjugate and 
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(pressure autospectra scaling functions from eq (5)), 
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If the pressure autospectrum does not vary significantly over the region of interest, W hh ((0) = W n ((0) = 1.0, 
and eq (1 1) simplifies to 


M a >) = 
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The term inside the bracket in eq (11) (or eq (12) for uniform pressure spectra) may be considered the squared 
transfer function between the total response at location j and the entire diffuse pressure field: 
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Note that H , is not an unchanging characteristic of the system since the scaling functions W depend upon 

the generally non-uniform patch autospectra, which change at different times in the launch trajectory. However, if 
the random pressure may be considered uniform over the surface of interest, H a , p is characteristic of the system 

independent of input pressure or the output acceleration. 


Eq (13) may be used to calculate 


Ajico) 


if the reference pressure 


ref 


and scaling functions 
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are known. 


Conversely, an effective reference pressure ma y be calculated if the acceleration PSD is specified and ^ 

is set to the identity matrix. The locations ^ selected as reference points on a skin without mounted components 
should be minimized, and chosen as far as possible from zone boundaries, joints, or local skin features such as 
cutouts, or doublers (which could distort the results). 


Now consider the case of a component mounted directly to the outer skin of a launch vehicle. The applied 
pressure does not change, but the transfer function in eq (13) developed for the bare skin must be generated for the 

component-loaded skin. Eqs (7) - (13) are applied at a response location using the modes of the component- 
loaded skin, and eq (13) becomes 


= h„,V w rVw h;,= 


A t (co) 


( 14 ) 
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where the tilde denotes the component-loaded acceleration/pressure transfer function and skin response. If the 
pressure is unknown, but the acceleration on the unloaded skin is known from measured flight data, the component- 
loaded response at any point t may be obtained by eliminating P f (co) from eqs (13) and (14): 


A f (cd) = Aj(oj)- 




(15) 


where the acceleration/single-patch pressure transfer functions H a lp ( CO ) in eq (13) for the bare skin are obtained 
from eq (7), and H a /p ( CO ) in eq (14) for the component-loaded skin from 


A',. M 


k=\ m= 1 


~ (J) film fikm 


co 1 -co 2 + i2A„ co„ co 

m ” m m 


(16) 


Equation (15) defines the RMM for predicting the response of a component-loaded skin structure when the bare 
skin response is known, and the acceleration/pressure transfer functions for both bare and component-loaded skins 
are available from modal analysis of the respective finite element models, or from modal testing. The frequency- 
dependent ratio of transfer functions in (15) replaces the Barrett scaling factor given by 

w 

B = n - — (17) 

w n +w c 


where W n and VV are the weights of the bare and component-loaded structures, respectively. 


Force and moment responses at specified interface elements (e.g., CBUSF1) may be obtained in the same fashion 
by replacing the first mode shape term in eq (16) with the modal forces and moments obtained in a NASTRAN 

RESTART in SOL 103. The ® term in the numerator of eq (16) is also dropped. The expression for the response 
force (or moment) at location ^ is similar to that for acceleration in (15): 

I ~ |2 


F cl (co) = Aj(co)- 



{co) 


|2 


(18) 


H ai/p {co) 


where F q (co) is the force (or moment) PSD at location q , and H f /p is the transfer function between the total 
force (or moment) at location q and the pressure: 




= H f VwrVwH t f - F ‘ ,ico) 
4 4 PM 


(19) 


The individual transfer functions ( CO) in H f between the force at location q and each single-patch 

pressure on any patch b are given by 
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( 20 ) 
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V. Validation 

Validation of the APTF method was provided in a previous technical presentation 14 . Therefore the RMM is 
featured in the validation section of this work. Examples will be supplied that overlay predictions from the direct 
APTF calculation with predictions from the RMM calculation. Furthermore, results calculated using the combined 
uncoupled transfer functions will also be compared to results calculated with the full-up coupled system transfer 
functions. These demonstration examples are currently in work. 

VI. Conclusion 

This paper presents an expansion upon successful validations of the APTF and RMM methods. The 
accomplishment provides a simple, but powerful, database analysis tool for the MSFC vibroacoustics team to 
estimate vibration responses at equipment mounting locations instead of using heritage processes that typically 
require running detailed analysis models and/or obtaining/processing substantial ground and flight test data. This 
puts cross-cutting capability into the hands of the propulsion and vehicle system departments to provide input 
vibration environment requirements for a new launch vehicle program. 
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The design and theoretical basis of a new database tool that quickly generates vibroa- 
coustic response estimates using a library of transfer functions (TFs) is discussed. During 
the early stages of a launch vehicle development program, these response estimates can 
be used to provide vibration environment specification to hardware vendors. The tool ac- 
cesses TFs from a database, combines the TFs, and multiplies these by input excitations 
to estimate vibration responses. The database is populated with two sets of uncoupled 
TFs; the first set representing vibration response of a bare panel, designated as H s , and 
the second set representing the response of the free-free component equipment by itself, 
designated as H c . For a particular configuration undergoing analysis, the appropriate H s 
and H c are selected and coupled to generate an integrated TF, designated as H s+C . This 
integrated TF is then used with the appropriate input excitations to estimate vibration 
responses. This simple yet powerful tool enables a user to estimate vibration responses 
without directly using finite element models, so long as suitable H s and H c sets are defined 
in the database libraries. The paper discusses the preparation of the database tool and 
provides the assumptions and methodologies necessary to combine H s and H c sets into an 
integrated H a+C . An experimental validation of the approach is also presented. 


Nomenclature 


Symbols 

A 

Acceleration power spectral density 

a 

Acceleration response 

f 

External force 

H 

Transfer function 

k 

acoustic wavenumber 

M 

Number of retained modes 

N p 

Number of pressure patches 

P 

Pressure spectral density 

P 

Pressure spectral density magnitude 

Pref 

Pressure reference spectrum 
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p 

Pressure 

r 

Distance from field point to response point 

W 

Pressure scaling functions 

a 

Areas associated with the degrees of freedom in a patch 

r 

Matrix of pressure field spatial function 

7 

Pressure field spatial function 

c 

Linear viscous damping ratio 

<t> 

Eigenvector 

LO 

Circular frequency 

Superscript 

c 

Component 

s 

Structure 

Subscript 

b 

Interface points 

e 

Excitation points 


Pressure patch indices 

l 

Degree of freedom index associated with single pressure patch 

m 

Mode number index 

N p 

Total number of pressure patches 

P 

Non-interface points 

£ 

Pressure patches 

Conventions 

0* 

Complex conjugate 

[] f 

Complex conjugate transpose 

Acronyms 

DAF 

Diffuse acoustic field 

LFP 

Large footprint 

MFP 

Medium footprint 

MSFC 

Marshall Space Flight Center 

NESC 

NASA Engineering and Safety Center 

RMM 

Response Matching Method 

RMS 

Root mean squared 

RPTF 

Response to Pressure Transfer Function Method 

PSD 

Power Spectral Density 

SFP 

Small footprint 

TF 

Transfer function 


I. Introduction 

P roviding suitable vibration requirements for launch vehicle equipment is a challenge faced by each new 
program. NASA’s Engineering and Safety Center (NESC) has identified methodology for producing vi- 
bration response estimates for mass-loaded vehicle panels as one of the top risks for new vehicle programs, 
observing that several different methodologies are currently in use with the potential to generate differing 
response predictions. Program managers need realistic estimates for vibroacoustic response environments so 
that vibration environment specifications can be developed for the design of hardware. Accurate prediction 
of vibration environments in early stages of a new program helps to assure adequacy of the hardware to with- 
stand the expected environments without specifying overly severe environments. Specifying environments 
that are overly conservative can increase financial and schedule costs as well as add weight to the vehicle 
design. Further, imposing less severe environments may result in hardware tests that can be performed more 
readily and with reduced cost. 

Generating accurate vibroacoustic predictions has historically required detailed analysis models and/or 
obtaining and processing substantial ground and flight test data. All of these activities require far more time 
and resources than are typically available at the start of a launch vehicle development program. Consequently, 
engineering organizations engaged in a new vehicle design need a vibration environment estimation process 
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that is quickly implemented and easily used. This paper describes a new tool that is being developed to 
fill this need. The methodology described here can quickly produce reliable vibration response estimates 
without requiring detailed analysis or test data. 

II. Background and Motivation 

Different organizations in the vibroacoustics community apply a variety of methods to estimate vibration 
environments for vehicle exterior panels. These environments are induced by fluctuating pressure fields on the 
panel from launch acoustics or flow- induced turbulent boundary layers during ascent. A review of the various 
methods in 2009, with concurrence by the NESC, indicated that when a panel is not loaded with mounted 
equipment, the different approaches returned results that were similar and consistent.^ Conversely, when 
vibration environments were considered for exterior panels with equipment attached, commonly referred 
to as mass-loaded panel systems, the predicted environments varied widely with little consistency across 
these same organizations. Such wide dissimilarity in the predicted environments throughout the community 
creates a high degree of uncertainty, and the NESC carried this finding forward as a risk for new launch 
vehicle (LV) development EE The risk has been addressed by the authors and othersPHSl 

In response to the NESC’s identification of a critical technological need, the authors have presented a 
pair of validated methodologies for calculating both vibration response and dynamic loads for equipment 
mounted to vehicle exterior panelsE! 7 The validation was based on acoustic response tests of a flight-like 
vehicle panel configured with different equipment masses. The first methodology is known as the Response to 
Pressure Transfer Function (RPTF) method. RPTF applies a correlated pressure field across the surface of 
a vehicle panel. In this application, a diffuse acoustic field (DAF) pressure assumption is adopted. Pressure 
spectra provided as RMS sound pressure levels are first converted to pressure autospectral density. Then 
cross-spectra associated with the pressure field excitation are calculated according to the best fit for a DAF. 
To implement this method, the excited vehicle panel surface must be subdivided into pressure patches. The 
distance between the center of each patch is selected based on the wave number required to calculate cross- 
spectra phasing between each patch. An acceleration response transfer function is first described in terms of 
an excitation on a single patch, but the transfer function is fully developed to include the spatial correlation 
of the pressure field across the many patches of the vehicle panel. 

The second method, the Response Matching Method (RMM), provides the basis for much of the theo- 
retical development presented here. The key difference between RMM and the present methodology is that 
previous versions of RMM have only been able to accept fully coupled TFs. The methodology presented in 
Section |IV| calculates the same response in terms of the uncoupled TFs of the bare vehicle panel and the 
attached component alone. This formulation enables panels and components to be combined in any desired 
configuration without the need for TFs specific to each configuration. 

RPTF and RMM have been developed such that they each accept as input one of the two most common 
forms of data available in the historical archive for past and present LVs: 

1. Zonal sound pressure levels are provided as standard requirements during the development phase of a 
new vehicle design and can be used as inputs to the RPTF method. 

2. Acceleration measurements on large-acreage exterior panels are often collected during the development 
flight tests of a new vehicle. Autospectral densities are often available either directly from a test flight 
or as an adjusted estimate from a similar vehicle design. These may be applied as inputs for predicting 
mass-loaded panel responses using RMM. 

Software that implements RPTF and RMM already exists, but requires that the analyst has access to 
the complete coupled form of the system equations for the mass-loaded panel of interest. The theoretical 
development in Section |IV| shows how RPTF and RMM can be extended to solve the coupled system equations 
while using the uncoupled system equations as inputs. 

Table [l] compares RPTF and RMM in their two forms. Each method has a form that accepts coupled TFs 
and a form that accepts uncoupled TFs. RPTF and RMM in both their forms offer the ability to accurately 
predict random responses; not only acceleration but also velocity, displacement, panel and component stress 
and strain, interface force, or any other linear output expressed as a set of modal eigenvectors. The new 
development of allowing RPTF and RMM to accept uncoupled TFs allows the user to quickly and accurately 
evaluate vibration environments for a wide variety of component configurations and interface designs without 
the need for costly and time-consuming updates to a large system FE model. 

3 oflHl 


American Institute of Aeronautics and Astronautics 


Table 1. Methods for populating the launch vehicle database tool library. 


Method 

Transfer 

Functions 

Modal 

Requirements 

Type of Input 

Available Output 

RPTF 

Coupled 

Loaded Structure 

Pressure Auto Spectral 
Density 

Any Response Auto and 
Cross Spectral Densities 

RMM 

Coupled 

Unloaded 
Structure and 
Loaded Structure 

Unloaded Skin 
Acceleration Auto 
Spectral Density 

Any Response Auto 
Spectral Density 

RPTF 

Uncoupled 

Unloaded 
Structure and 
Component 

Pressure Auto Spectral 
Density 

Any Response Auto and 
Cross Spectral Densities 

RMM 

Uncoupled 

Unloaded 
Structure and 
Component 

Unloaded Skin 
Acceleration Auto 
Spectral Density 

Any Response Auto 
Spectral Density 


III. Description of Database Tool 


Extending RPTF and RMM to accept uncoupled TFs as inputs enables the development of a new 
database tool to predict the mass-loaded panel response for a new LV. The database tool allows an engineer 
who may not be well- versed in finite element methods to quickly visualize the narrow band vibration response 
of a given primary structure panel with mounted components selected from a library. Each component is 
integrated into the vehicle panel automatically, with no requirement to create a new FE model. The library 
of available components that can be integrated into the vehicle panel will span a realistic range of weight, 
centroud offset, and interface footprint. 

A typical example of a component library for a specific LV zone is shown in Table [2j The size of 
the footprint for each component is designated as small (SFP), medium (MFP) or large (LFP). When a 
component is selected, the database tool generates the coupling equations derived from the unloaded panel 
and component natural frequencies and modes and then applies the panel pressure or acceleration input to 
predict the mass-loaded panel response. If more detailed modal data for the component is known from either 
test or FE analysis, this information may be used instead of the pre-loaded library components. 


Prepare Bare 
External Vehicle 
Panel FEMS 


Use FEM based tools to construct two sets of uncoupled Transfer 
Functions and Load Database Tool with Zonal Sound Pressure Levels 





Panel Transfer 
Function H s (w) Set 

Interface Response/ 
Correlated 
External Pressure 
TF Set 


1 




Production Team uses data base 
tool to assemble transfer 
functions select Excitation 
Environments and produce best 
in class Response Estimates 


Assemble 

Bare External Vehicle 
Panel Transfer 
Functions H s (w) 




Prepare Equipment FEMS 


Equipment Transfer 
Function H c (ui) Set 


Component Response/ 
Interface Force 
TF Set 


Select Response Locations 
for Equipment Interface 
locations 


Select Frequency Resolution 
and Damping Assumption 
for Response Solution 


produce Frequency 


Assemble 
Component 
Equipment Transfer 
Functions H c (u>) 


Figure 1. Process flow for preparing the database tool 
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Table 2. Example of possible choices for components in the database tool library 


Component 

Designation 

Number of 
Interfaces 

Interface 

DOF 

Weight 

(lb) 

Centroid 
Offset (in) 

l-SFP-15-8 

4 

24 

15.00 

8.00 

2-SFP-30-8 

4 

24 

30.00 

8.00 

3-SFP-45-8 

4 

24 

45.00 

8.00 

4-SFP-15-11 

4 

24 

15.00 

11.00 

5-SFP-30-11 

4 

24 

30.00 

11.00 

6-MFP-30-8 

8 

48 

30.00 

8.00 

7-MFP-45-8 

8 

48 

45.00 

8.00 

8-MFP-60-8 

8 

48 

60.00 

8.00 

9-MFP-15-11 

8 

48 

15.00 

11.00 

10-MFP-45-11 

8 

48 

45.00 

11.00 

ll-LFP-45-11 

16 

96 

45.00 

11.00 

12-LFP-60-11 

16 

96 

60.00 

11.00 


Figure [l] shows the process flow for preparing the database tool library. The curved panel in this il- 
lustration depicts the 16 interface locations for a LFP component. The tool generates the coupled system 
response from uncoupled panel and component transfer functions derived from the pre-computed or mea- 
sured modes. The database tool preparation team generates the necessary and consistent data for the library 
from finite element models or measurements and defines default values for several user-accessible database 
tool parameters. 

The database tool is designed to anticipate the need to expand in scope and complexity as time passes and 
a LV program matures. In the beginning of a LV development program the primary structure designs may 
be known in rigorous detail. The component design and integration details may not be as well-defined. A 
set of default component options should be sufficient to quickly develop reasonable vibration environments, 
qualification test specifications, and preliminary cost estimates early in the program. As more detailed 
information becomes available, higher-fidelity FE models of the actual component and integration details 
are easily installed in the database for more precise and accurate estimates of the vibration environments. 

The database tool end-user does not create a model or perform an FE analysis. Instead, he or she makes 
reasonable selections from panel structures and components loaded in the library. He or she also selects the 
flight event of interest (e.g., lift-off, ascent). The selected flight event and the vehicle zone containing the 
selected panel will correspond to two sets of source environments from which the user may choose one: either 
the sound pressure level or unloaded panel acceleration spectral density. 

After making these selections and accepting the default setting for a damping schedule the database tool 
will calculate the coupled system transfer functions from the uncoupled primary structure and component 
modal characteristics. The theoretical basis for the coupling methodology is provided in the next section. 

The benefit of assembling the loaded panel transfer functions from the separate unloaded panel and com- 
ponent transfer functions is the ability to predict the loaded panel vibration environment without requiring 
an update of the full system modal solution. Thus, the computationally expensive modal solution of the 
vehicle is required only once for a vehicle design cycle, and components may be selected from the default 
database tool library or provided anew as the component design matures. But lunch is never free. To achieve 
reasonable accuracy in the final response prediction from decoupled transfer functions, more modes need to 
be retained for both the unloaded structure and the component models. (Preliminary data suggest that 



American Institute of Aeronautics and Astronautics 


the mode retention requirements are especially large for the component models.) Further refinement of the 
database tool will incorporate well-known component mode synthesis techniques to alleviate some of the 
added computational burden. 


IV. Methodology 


We begin by deriving an expression for the response of a structural panel with an attached component. 
The relation is written in terms of the individual transfer functions associated with the panel and compo- 
nent. Since the database tool contains only individual panel and component transfer functions, writing the 
combined transfer function in this manner is essential. Similar derivations have been presented by Majed, et 
al. 6 Next, the pressure field acting on the structure is considered. By conveniently expressing the pressure 
field in terms of known quantities, the acceleration spectral density of the component loaded panel can be 
found. It is this expression, Eq. (20), that underpins the core computations of the database tool. 


IV. A. Response of Structural Panel with an Attached Component 

Consider a structural panel with an attached component as shown in figure [2| The panel is subject to force 
excitation at e locations and has interface points with the attached component at b locations. A frequency 
dependent transfer function, R(w), is defined simply as the ratio of the acceleration response to the input 
(e.g., H = a/F ). In this case, the input is an external force; pressure inputs will be considered subsequently. 
Given the definition of the transfer function, the response of a structure is generally given by 


a =[H]{{}, 


Cl) 


where the explicit frequency dependance has been dropped for convenience. Partitioning Eq. |T] to consider 
the response at b interface locations and p non-interface points gives 




(2) 


From a free-body diagram of the component, the forces at the interface are related to the accelerations by 

f b = -[R b c h ]- 1 a b . (3) 

Combining Eqs. <[2j) and <[3j) and solving for a b we may express the acceleration at the interface in terms of 
the interface transfer functions and the external force 


a b 


Wi 


bbi 


[Hi 


bbl 


[H&r'W Jf e 


( 4 ) 


Now consider points p that are located on the panel, but not at interface points. The acceleration response 
of such points is given by 


■=[H 1 


-pej 


\P-pb\ *b- 


Using Eqs. ([3]) and ^ as well as acceleration continuity, the response at points p 


a„ = 


[H s pe ] - [H; b \ [[!%] + [H^]}- 1 [H s be 


u. 


( 5 ) 

( 6 ) 


Eqs. ( |4| ) and. @ will be used in Section [IV.C| to find the response of the component loaded panel subject to 
random pressure excitation. 

Now consider the response at a single point p due to a single force excitation e. For linear structural 
systems, the complex scalar elements of the force/acceleration TF matrices shown in Eqs. 0 and 0 can 
be expressed as a summation of natural frequencies and modes 


M 

Hpe = 

m = 1 


(<4 


^ 0pm 0em 

- w 2 + ’ 


( 7 ) 
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Figure 2. Schematic of structural panel (s) with attached component (c) 


where the TF elements given by Eq. (J 7 | take the same form for responses at interface and non-interface 
points. This form also applies regardless of whether one is assembling matrices of bare structure TFs or 
component TFs (provided, of course, that the appropriate structure or component natural frequencies and 
modes are used). 


IV. B. Pressure Excitation of Panel 

A random pressure held on a launch vehicle skin surface may be approximated by dividing the surface into 
patches, or regions of uniform pressure with no phase offsets within the patch. The size of each patch must be 
chosen to be small enough to justify the assumption of uniform (though dynamically varying) pressure with 
zero phase offset across the patch. The selection of the appropriate patch and element size is not included 
in this monograph, but is addressed in a separate studyP 

The pressure may be defined as a stationary Gaussian random field with spatially varying autospectral 
density. The pressures on any pair of patches may be correlated, exhibiting a non-zero cross-spectral density 
between them. The random pressure field is thus a Hermitian matrix of spectral densities of dimension N p , 
the total number of pressure patches. The pressure autospectra occur on the diagonal of the matrix. The 
cross-spectra appear off of the diagonal. The random pressure field on all patches may be written as 



Pll 

P12 ■ 

PlN p 

P(u) = 

P21 

P22 

■ P 2 N,, 


. p Npl 

Pn p 2 ' 

■ Pn p n p 


( 8 ) 


where Pjk = -Pjb. If spatial functions 7 (w,r) are defined to relate the autospectra to the cross-spectra, 
Eq. (|8j) may be written as 


7 hAi 

712 A2 

llN p P\N p 

721 P12 

722-P22 

r )2N P P2N p 

'7n p iPin p 

lN p 2P2N p ■ 

■■ 7Af p Af p -PAf p Afp 


where Pjk = yjPjjPkk- As will be discussed shortly, in the case of a DAF, 7^(07 r^) = 1 along the 
diagonals; however, the 7 terms have been included in the diagonal terms in Eq. © for generalization. The 
expression for Pjk satisfies an inequality requirement on the coherence which states that 


0 < 


\ p jk(u )\ 2 

Pjj ) Pkk ) 


< 1.0. 


(10) 
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For a DAF, the spatial functions may be expressed as 


7 jfcM r jk ) 


s in ( hr jk) 
krjk 


( 11 ) 


where the patch centroids are constrained to lie on the curved skin surfaces. When j = k the spatial functions 
coincide with the patch autospectra and the gamma approaches unity by L’Hopital’s Rule. Finally, the 
components of the pressure matrix may be expressed as products of frequency-dependent scaling functions, 
Wjk and an arbitrary reference autospectrum pressure, P re f, he., 


— ^ijk ll j/c Pef ) 


(12) 


Substituting Eq. ( 12 ) into Eq. (|9| 
lift-off environments in this study 


we obtain an expression for the diffuse-field pressure model used for the 


or more compactly, 


where for a DAF, 
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(13) 


(14) 


(15) 


It is important to recognize that all of the pressure auto and cross-spectra discussed are related to the 
known reference spectrum, P re f through the scaling functions. If the acoustic field is not diffuse and the 
scaling functions are unknown, then the skin structure under consideration must be limited to regions with 
similar pressure autospectra and must be set to unity. The scaling functions may also be grouped into fewer 
distinct sets for large regions showing uniform pressure autospectra. That is, many of the functions may 
be assigned the same value for various zones of the launch vehicle even if the zone is subdivided into many 
patches. 


IV. C. Mass Loaded Panel System Response Due to Pressure Field Excitation 

The acceleration response of a panel at non-interface points p due to a vector of patch pressures, p e is 


= [H;t c ] { P j. 


(16) 


where IP+ C denotes the acceleration/pressure TF of the component loaded panel. From Eq. |6j), it follows 
that Hpf c can be written in terms of acceleration/pressure TFs for the bare structure and component as 


[H s p t c ] = [H; e \ - [H; b \ M + [H^r 1 [HZ S ] . (17) 

Acceleration/pressure TFs can be related to acceleration/force TFs by considering the pressure on a patch 
as an array of distributed point forces. To then convert an acceleration/force TF to an acceleration/pressure 
TF it is only necessary to multiply by the surface area associated with each of the point forces on a patch. 
(Computing the surface area associated with a given node is a calculation that can be readily performed 
with a finite element software.) Using this TF conversion and Eq. ([?|, it is possible to write the individual 
elements of [H pe \ corresponding to response point p and patch e 


M 


H„ 


= £ 


, -UJ 2 {4>l m } {Qt} 

pm (M - w 2 + 2 i? m w m w) ’ 


( 18 ) 
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where {<f>im} is a partial eigenvector containing the elements corresponding to the degrees of freedom asso- 
ciated with a single patch e and {cq} is the vector of areas associated with those degrees of freedom. 

Turning attention to the acceleration PSD response of a panel, the total response at a single point 
p includes the autospectra from the pressures on all of the patches and also from non-zero cross-spectra 
between any two patches,® or 

N p N p 

A P = IlIl H *pe j H pSk P jk . (19) 

j k 


Using Eqs. (14) and (19), the PSD response of a component loaded panel at point p can by expressed in 
matrix-vector form for computational efficiency as 


A ; +c = { H a P i c } m* tn [wY{H;t c yPre S , 


( 20 ) 

where is a lxN p vector and N p is the number of pressure patches. The individual elements of 

{H p + C } are given by Eq. (18). Similarly, the PSD response of the bare structure is 


a; = {h; £ } \w]* [r] \w\^{H; e } ] p ref . 


(21) 


If the pressure field applied to a component loaded structure is unknown, but the acceleration PSD response 
of the bare structure to the same pressure field is known, it is possible to combine Eqs. (20) and (21) to 
obtain a relation that is independent of pressure 


AS+C _ AS 

A P ~ ^ P 


{H;+c}[wf>[T} \W)Hh£ c Y 

{H° )e }\wY [r] [wY{H» )e }' 


(22) 


V. Equivalence of the RMM and the RPTF Method 

As discussed in Section [TTJ NASA’s historical archive of source environments data for past and present 
launch vehicles contains data as either zonal sound pressure levels or acceleration autospectral densities 
corresponding to specific locations on the launch vehicle surface. The database tool will enable the choice 
of either of these inputs. When sound pressure levels are used, the mass-loaded panel environments will be 
calculated using RPTF. If the source environments are in the form of acceleration autospectral densities, then 
the output will be calculated using RMM. It is important to note that these are mathematically equivalent 
methodologies; their only substantive difference is the type of input that they accept. Figure [3] shows the 
equivalency of the RPTF method and RMM. Figure [3] (a) shows the source environment as a pressure power 
spectral density. This pressure environment was propagated through a bare panel using RPTF and the 
resulting acceleration power spectral density is shown in Figure [3] (b). Figure [3] (c) shows the acceleration 
power spectral density of the mass-loaded panel calculated using both methods. The RPTF solution (solid 
line) started with the pressure spectral density in Figure [3] (a) and used the uncoupled transfer function 
formulation of the method to calculate an acceleration power spectral density for the mass-loaded panel. 
The RMM solution (dashed line) started with the bare panel acceleration power spectral density curve 
shown in Figure[3](b) and also used the uncoupled TF formulation to calculate a mass-loaded panel response 
curve. As expected, the two curves are identical. All of the presented source environments and responses 
correspond to the direction normal to the panel; however, the equivalency of the two methods in the other 
two axes can also be shown. 

While the mathematical equivalency of the RMM and RPTF methods has been shown, it is important 
to note that in practice, the two methods often will not predict equivalent responses. This is due to the 
manner in which source environments are created and stored in the historical archives. Pressure spectral 
density source environments for lift-off acoustics are inherently smooth and broadband in nature. They 
are also relatively repeatable from launch to launch. Thus, when these environments are enveloped and 
stored in archives, only a modest amount of conservatism must be added to the spectrum to account for 
the uncertainty in the environment. This is not the case for acceleration power spectral density source 
environments. Due to the natural modes of the vehicle panels, these source environments tend to have many 
more discrete peaks. Furthermore, the frequency at which these discrete peaks occur will have some flight- 
to-flight variability. In the process of enveloping these environments for storage in the archive, a relatively 
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Figure 3. Equivalence of the RPTF method and RMM. (a) Pressure power spectral density source environment, (b) 
Equivalent source environment in terms of acceleration power spectral density, (c) Acceleration power spectral density 
of the mass-loaded panel calculated using both methods. 


large amount of uncertainty is applied. Therefore, in practice, the use of RMM and acceleration power 
spectral density source environments will typically result in higher response spectra than will result from 
using RPTF and a corresponding pressure spectral density source environment. While both types of source 
environments may not always be available for a given launch vehicle panel of interest, if the analyst has a 
choice, it will generally be more conservative to use RMM with acceleration power spectral density source 
environments. This more conservative approach is typically warranted early in a launch vehicle development 
program. As a launch vehicle program matures, it may be possible to remove some conservatism from 
component requirement specifications by using RPTF with the appropriate pressure spectral density source 
environments. 


VI. Implementation and Validation of Uncoupled TF Formulation 


To prove the viability of the database tool concept it is necessary to show that the response of a mass- 
loaded panel generated by the uncoupled TF formulation is approximately equivalent to the response gen- 
erated by the coupled TF formulation. To show this approximate equivalence a finite element model of 
a flight-like panel with and without attached equipment assemblies was developed. A modal analysis was 
performed on these models using MSC/NASTRAN and the corresponding TFs (i.e. , H s and H s+C ) were 
calculated using the appropriate equations (either Eq. ([7]) or Eq. (18). A free- free modal analysis of the 
attached equipment assembly model by itself was also performed to calculate the necessary component TFs 
(i.e., H c ). The finite element models were designed simulate a configuration recently ground tested in the 
acoustic chamber at MSFC. A companion paper by the authors details this test series and presents the 
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Figure 4. Oblique view of FE model of panel and surrounding support structure. 


currently available data? 9 Figure [4] shows the finite element model of the panel and the attached equipment 
assemblies. The surrounding wall structure was also modeled and is shown here. This model contains 60,841 
nodes and 57,313 elements. Checks were made to ensure the natural frequencies and mode shapes predicted 
by the models showed good agreement with those found in testP 

The test panel was subjected to acoustic noise excitation energies closely simulating the liftoff event and 
roughly approximating ascent flight events. This acoustic pressure spectrum was then used as an input to 
the RPTF method. The coupled TF and uncoupled TF formulations of the RPTF method were then used 
to calculate the response of the mass-loaded panel at an interface point. The results are overlaid in Figure[5] 
Note the excellent agreement between the responses across the entire frequency range. The agreement is not 
perfect, however. It is suspected that this discrepancy is due to the approximate nature of the uncoupled 
formulation used here. Other component mode synthesis approaches (e.g., the Craig-Brampton method 10 ) 
would likely result in a more exact match. Implementing such approaches will be the subject of future 
work, but the agreement evidenced by figure [5] is more than adequate for the purpose of quickly developing 
vibroacoustic response estimates early in a vehicle program. 

Responses were also calculated for a non-interface point on the panel. A point on the model was chosen 
to correspond to a location at which test data were available. The acceleration response spectra from the 
uncoupled RPTF formulation, the coupled RPTF formulation, and the test data are overlaid in figure [6] 
The coupled and uncoupled formulations show excellent agreement with each other and the test data. This 
serves as further validation of the coupled RPTF formulation and the first ever validation of the uncoupled 
RPTF formulation. Since the uncoupled RMM formulation was shown in Section [V] to be mathematically 
equivalent to the uncoupled RPTF formulation, the validation of the uncoupled RPTF formulation implicitly 
validates the uncoupled RMM formulation as well. Note that the damping schedule used to calculate both 
analytical responses is a damping schedule identified from the test data using a software tool developed at 
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Frequency (Hz) 


Figure 5. The response of a mass-loaded panel at an interface point calculated using the coupled TF and uncoupled 
TF formulations of the RPTF method. 


MSFC known as DampID. More information concerning DcimpID is presented in Smith, et al. 9 

VII. Conclusions 

This paper presents and validates the uncoupled TF formulations of the RPTF and RMM methods. This 
enables development of a simple, but powerful, database analysis tool for the MSFC vibroacoustics team to 
estimate vibration responses at equipment mounting locations. With this database tool, it will no longer be 
necessary to use heritage processes that typically require the development of detailed analysis models and/or 
the collection and processing of substantial ground and flight test data. This puts powerful capability 
into the hands of the propulsion and vehicle system departments to provide input vibration environment 
requirements for a new launch vehicle program. In addition to improving the uncoupled formulations with the 
use of well-known component mode synthesis techniques, future work may involve expanding the methods to 
accept other types of source environments as input. For instance, the aero-fluctuating pressure environment 
associated with vehicle ascent is often modeled with what is known as a Corcos model Implementing such 
a model into the database tool for use with the RPTF method would provide even greater capability for 
predicting environments across all flight regimes. 
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Figure 6. Acceleration response spectra from uncoupled RPTF formulation and the coupled RPTF formulation com- 
pared with test data. 
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